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We explore fermion RG block-spinning transformations on the lattice with the aim of studying 
the IR structure of gauge theories and, in particular, the existence of IR fixed points for varying 
fermion content. In the case of light fermions the main concern and difficulty is ensuring locality 
of any adopted blocking scheme. We discuss the problem of constructing a local blocked fermion 
action in the background of arbitrary gauge fields. We then discuss the carrying out of accompa- 
nying gauge field blocking. In the presence of the blocked fermions implementation of MCRG is 
not straightforward. By adopting judicious approximations we arrive at an easily implementable 
approximate RG recursion scheme that allows quick, inexpensive estimates of the location of con- 
formal windows for various groups and fermion representations. We apply this scheme to locate 
the conformal windows in the case of SU(2) and SU(3) gauge groups. Some of the reasons for 
the apparent efficacy of this and similar decimation schemes are discussed. 
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1. Introduction 

The dependence of the phase structure of non-Abelian gauge theories on their fermion content 
is a question of obvious field theoretic and physical interest. In view of anticipated discoveries at 
the LHC this question has recently come to the fore with some urgency. For sufficiently small num- 
ber of fermions there is UV asymptotic freedom (AF) and IR confinement. It is also known that for 
sufficiently large number of fermions asymptotic freedom and confinement are lost. It is generally 
believed that there is an intermediate range with an UV AF fixed point but no confinement; instead 
there should be an IR fixed point resulting in conformal IR dynamics. The dependence of these 
'conformal windows' on the gauge group, the number of fermion multiplets and their represen- 
tation, as well as the structure of the associated IR fixed points, i.e. their relevant directions and 
anomalous dimensions, are questions of central importance for informed model building "beyond 
the standard model". For example, the viability of current walking technicolor models depends on 
such detailed knowledge. 

Several approaches have been employed over the last few years in order to investigate these 
issues (see [1] - [2] for review and references). One popular method is to adopt some definition of a 
running coupling, e.g. the 'Schroedinger functional' definition. Since any definition of a 'coupling' 
outside the perturbative regime is necessarily arbitrary, any such choice is scheme-dependent. Cau- 
tion should then be exercised in interpreting results. Another approach is to compute the spectrum 
and/or some judiciously chosen physical observable and track their deformation as a function of 
the number of fermions and their masses. Such spectral studies should, in principle, yield a clear 
outcome, but they are very expensive to carry out in practice to the degree of accuracy required for 
unambiguous results. Furthermore, it is not always clear what the expected behavior should be. 

A direct method for getting at the phase structure would be implementation of the Wilsonian 
Renormalization Group (RG). Successive RG blocking transformations determine the effective 
action RG flow and its fixed points. Conceptually, this is the most straightforward and unambiguous 
approach. Various RG-based studies have been carried out recently for the most part using the two- 
lattice matching MCRG method [4] - [7]. There are, however, special challenges that must be 
met when attempting to implement RG blockings and MCRG methods in the presence of light 
fermions. 

In the following we will examine some of the issues involved in performing RG transforma- 
tions with fermions. We will examine the question of preserving locality of the resulting action 
after an RG blocking step of light fermion degrees of freedom. For full RG implementation such a 
step must of course be accompanied by a gauge field blocking step. Again, how to do this in prac- 
tice in the presence of fermions, in particular how to implement MCRG, is not straightforward. To 
get some insight we will adopt some judicious approximations that lead to an easily implementable 
approximate RG recursion scheme [8]. This scheme, which proves rather effective, is then applied 
to locating the conformal windows in the case of the SU (2) and SU (3) groups. Some or the reasons 
for its apparent effectiveness are discussed in the last section. 

2. Fermion RG blocking 

In devising RG blocking transformations for light fermions the main concern is maintaining 
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reasonable locality of the blocked action after each blocking step. This is a question of interest in 
it's own right apart from issues of practical implementation. For free fermions a blocked action 
can be defined along the lines below (eq. (2.3)-(2.4)) that can be shown to maintain locality [9]. 
This has been used to argue [10] the locality of determinant rooting for free fermions. This locality 
demonstration, however, cannot be extended in the presence of an arbitrary gauge field background. 
To begin exploring this issue consider the Wilson operator in background gauge field U : 

Do(U) = \-kM[U], (2.1) 

where the only non- vanishing elements of matrix M are: 

M n{n+fl) [U] = (l-r fl )U tl (n), M n[n _ (l) [U} = (l + r ll )Ul(n-fL). (2.2) 

A fermion RG blocking step specifies a transition from the fermion degrees of freedom on the 
original lattice A to 'thinned-out' degrees of freedom on a lattice A'. With A a lattice of spacing 
a, let AW denote the lattice of spacing la. We take A' = A' 1 ) US, where S is a subset of A to be 
specified later. 

We introduce variables on A' = A^ 1 ' U S by 

J Dij/Dy exp(— xj/Doy/) 

= j D\j/DYD^ ([) DW [[ Uxp{-\j/DoY-a{^ {[) -YQ^) {^ {[) -QV)} ■ (2-3) 

The rectangular matrices Q : A — > A' are given by 

d d 
QnmYm = Yn + C L M n («+/t) P] Y(n+fi) +^1^ („_£) [U ] Y{n-p.) (2-4) 

= {$nm + £Qnm) Ym ■ (2-5) 

Note that a — > °° gives a 5-function definition of the *F (1) 's, i.e. = (QY)n- Integrating out the 
original variables one has on the new lattice A': 

/ DYD\ff exp(— vrDovO 

J A 

= Deter 1 J D»F (1) D»F (1) exp{-a^ (1) »P (1) + a 2 ^ (1) eGe t ^ (1) } (2.6) 

with 

G l = [D Q + aQ?Q] . (2.7) 

In general, with light fermions, G is a very non-local propagator. 

We now choose the set S and the decimation parameter £ in (2.4) so that in the action 

-xj/Doy/- a (^ 0) - YQ^) (^ (n - Qy) 

in (2.3) one cancels the nontrivial part of Do: 

Y[KM-a^(Q + &)] ¥ = 0. (2.8) 
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(a) (b) 

Figure 1: (a) The set S, as a subset of A, is the set of sites reached from the interior site of each 2x2 
hypercube by displacement by a in each lattice direction; (b) fermions at sites in S interact with fermions on 
AW and among themselves by 'fat links' consisting of all possible 'staples' of length 2a, a subset of which 
is indicated here by broken lines. 



This is accomplished by setting £ = ff/a, and taking the set S as shown in Fig. 1(a), i.e. S, viewed 
as a subset of A, is the set of the interior sites of the boundary plaquettes of the 2x2-hypercubes 
that are the elementary d-cells of A^ 1 '. The propagator G for the the original variables \jf in (2.3) is 
then given by 

G-' = 1(a) + t-ffQ > . 1(a) = ( (' + "> 5 "'" \ A ' , (2.9, 

a o nm n,m £A\A 

and the action in (2.6) from the integration over the original \j/ becomes 

a 



1 + a 



K QGQ 



W>. (2.10) 



(2.10) is manifestly local on A' = US. In particular, note that from (2.9), taking a large, 
interactions induced by G may be expanded in a hopping expansion in ^ and resummed. In fact, 
in the limit a — > °°, i.e. for the 5-function definition of the blocked variables (cf. above), G becomes 
ultralocal, and in fact [/-independent, and the action goes to 

[l - K 2 Q&] . (2.11) 

After this blocking step we end up with fermions on the boundary of the blocked cube, i.e. a 
form of exact 'potential moving' fermion blocking albeit on a non-hypercubic lattice. The fermion 
action is still bilinear in the blocked fermion fields which interact via 'fat' gauge field bonds (Fig. 
1(b)). The gauge field bond variables are yet to be blocked. 

It is important to note that if one starts integrating out gauge field bonds inside each blocked 
hypercube, four-fermion and higher-fermion interactions will in general be induced (but of course 
at the scale of the blocked cube). Ultralocality of G in this connection completely avoids potential 
difficulties from the DetG~' factor in (2.6) in performing any such accompanying gauge field 
blocking. 

We have devised a fermion blocking step in a general gauge field background that results in a 
manifestly local action. How to perform successive fermion blocking steps in a fixed gauge field 
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background under the requirement that the action be manifestly local at the scale of each step is 
an interesting question in its own right. The following is a sufficient condition: for the n-th step 
implement the above blocking procedure integrating out the fermions on the sites satisfying the 
conditions: 

X[ = 2"~ 1 ra ; - , nti 6 Z, i = 1, . . . ,d , 

£jc,- = 2 n l + 2 n - 1 , leZ. (2.12) 

i 

Then the n-th step blocked fermions interact again via bilinear interactions through gauge field fat 
links of maximum length 2" (in units of a). One can actually do better if one allows maximum 
length 2" +1 and performs a two-stage blocking of fermions interior to each blocked hypercube. 
The general lesson, though, is that one cannot achieve a manifestly local blocked fermion action 
residing on a purely hypercubic lattice, i.e. without some non-empty analog of the set S above (cf. 
Fig. 1). This seems to accord with previous studies of fermion actions [1 1]. 

3. Full RG blocking 

For the physical questions we are here interested in, i.e. determination of the RG flow and its 
fixed points, RG blockings of both fermions and gauge fields must of course be implemented. A 
local blocking step for fermions such as above must then be accompanied by a gauge field block- 
ing step. It is of course possible to define a variety of gauge field blocking transformations, and 
several instances have been employed the literature, e.g. Swendsen, DSB or HYP-type transfor- 
mations. Having defined the blocked gauge field variables, however, actually implementing the 
gauge blocking can generally be done only numerically. In the case of pure gauge theories, MCRG 
provides a tried method. Generally, given a starting action one generates a set of configurations by 
MC which are blocked. An effective action on the blocked lattice is then determined by measuring 
its coupling on the blocked ensemble. Effective actions that reproduce expectations of observables 
over a wide range of scales can in principle be so constructed [12] - [13]. This procedures may be 
repeated over several blocking steps. If one is not interested in the effective action flow but only in 
observables, the two-lattice matching method [5] may be used. Here one generates two ensembles 
of configurations, starting from two different couplings in the original action; a matching procedure 
is then applied to them by measuring observables on them. 

In the presence of fermions, however, it is not clear exactly how such MCRG procedures are to 
be applied. One cannot generate a set of configurations to be blocked without first integrating out 
the fermions completely. One may then apply MCRG to this ensemble. In particular the two-lattice 
matching method has been applied in this manner in [4] - [7]. Eliminating the fermions of course 
generates a very non-local action for the gauge fields due to the resulting fermion determinant in 
the gauge field background, and naive universality notions do not necessarily hold. This raises 
some issues. It is not immediately apparent that measurements on repeatedly blocked configura- 
tions generated from this non-local action are equivalent to measurements using the local effective 
action in terms of the blocked fermion and gauge variables resulting from local blocking steps. In 
particular, it is not clear that for fermionic observables naively written in terms of simple fermion 
degrees of freedom on the blocked lattice sites some kind of universality guarantees equivalence to 
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(a) (b) 



Figure 2: Approximations: (a) Merge the fermions on S from the exact blocking with the fermions on the 
hypercube edges while appropriately 'renormalizing' the hopping parameter of the latter, cf. text; (b) stan- 
dard potential moving approximation for gauge field blocking: move all interior plaquettes symmetrically to 
the boundary with appropriate weights. 



measurements with the local blocked action and observables in terms of blocked fermion degrees 
of freedom. Devising ways to test assumed equivalence of procedures, or designing new MCRG 
procedures for full fermion-plus-gauge field RG blocking still remain a largely unexplored area. 

Here, in order to get some insight, we opt for the modest goal of devising and trying out an 
approximate decimation scheme [8] that can easily be explicitly earned out. 

4. Approximate RG scheme 

To obtain such a scheme we proceed as follows. Eliminate first the fermions on S and com- 
pensate by adjusting the hopping parameter of the remaining fermions residing on the hypercubic 
lattice - this may be viewed as a typical fermion potential moving approximation to get a purely 
hypercubic action with 'renormalized' gauge interactions (Fig. 2). The renormalized hopping pa- 
rameter is determined by requiring that some long-distance gauge-invariant fermion correlation 
function remains invariant under the fermion blocking step. Thus, taking the 2-point function, 
\j/ n U[C nm ]Ym, (other choices, such as a meson-meson correlator, will do equally well) this require- 
ment translates into the exact relation: 

K(2a) = D(K(a))K(a) 2 , (4.1) 

where 

£ ( %(1 - Yn) (! - Yv) V [C nm ] Up (n)U v (n + A) Vm) 

D(K(a)) = ^— — , . (4.2) 

2, ( ¥n( 1 - Yn) U [C nm ] U^(n)U^(n + jl) y m j 

(4.1)-(4.2) provide a recursion relation relating the hopping parameters at lattice spacings a and 2a, 
which may, in particular, be solved to obtain the critical hopping parameter. One may then proceed 
to carry out the gauge field blocking step by whatever method is chosen. 

The above approximation maintains bilinear fermionic interactions at each successive blocking 
step. This suggests the following considerations. Assume one implements some RG blocking 
scheme. Denoting by U^ n \ \jf("\ the blocked variables at the n-th step, then, quire generally, 
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one has 

/ DU^Dy^Dy^ exp(-S { " ] - 
Jam ' 

= f DU^Df {n+1) Dxi/ in+1) exp(-4 +1) -4 ,+1) )- (4-3) 

If, furthermore, as after the above approximation, Sp and Sp +1 ^ maintain the same bilinear form 
on AW and A^" +1 \ respectively, then 

exp(-4 +1) )= / DU^F(U {n \u {n+ V) 
Jam 

with F(t/W,t/( B+1 )) specifying the blocking relation between and LA" +1 )'s. The ratio of 
determinants in (4.4) can now be expected to be reasonably local. Indeed, major part of non- 
locality is manifestly cancelled in the ratio. This can be seen from the formal expansion of the 
determinants in terms of gauge field loops: 

Detsf = n det (! - K |C| M[C]) , (4.5) 
c 

where M[C] denotes the path ordered product along the loop C of the matrices M(U) defining 
the fermion bilinear operator in Sp, the product in (4.5) being over all loops. All loops in the 
expansion of lnDet4" +1) then are contained in that of DetSp . The main contribution to the ratio 
comes in fact from loops inside each blocked hypercube. This suggest that one may approximate 
it by just such a set. 

4.1 Simplest RG implementation model 

These considerations lead us to an approximate implementation of RG decimations of the full 
system of fermions and gauge fields given by the following recipe. 

After a fermion blocking step, apply the approximation depicted in Fig. 2(a) combined with 
(4.1)-(4.2). Then the gauge field action on the spacing 2a lattice is given by (4.4). As just explained, 
after the cancellation of the denominator by the numerator in the determinant ratio expressed in 
terms of gauge field loops, we can approximate the remainder by the contributions of loops interior 
to each blocked hypercube - in particular, by just the single plaquette loop weighted by a parameter 
r/. (One may of course consider approximations involving more interior loops, e.g. 2x2 loops, 
and more decimation parameters.) It remains to specify and perform a gauge field blocking over 
each hypercube. This we do by standard potential moving decimation: every interior plaquette is 
moved to the hypercube boundary with the appropriate weights 2 cl ~ 2 (Fig. 2(b)). Inserting these 
steps in (4.4) one obtains 

exp(-4" +1) [[/[p(" +1 )]])=/ J] dU b <ta(l-**M\pV>]f* J] (exp(-4 n) ^[P (n) ]])) 2d " 

(4.6) 

giving the plaquette action for the gauge field ,Sg +1 ^[£/[//' !+1 )]] on lattice A(" +1 ) in terms of the 
theory on A^. In (4.6) G signifies the plaquettes pW in A*") tiling plaquette p( n+1 ) 



DetS 



(/<) 



DetS 



n+l) 



exp(-S, 



(4.4) 
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N L 


Mj 


SU(2) fund 


9 


11 


SU(2) adjoint 


2 


3 


SU(3) fund 


11 - 12 


17 


SU(3) sextet 


2 


4 



Table 1: Conformal window lower and upper bounds as obtained from the blocking recursion relations 
described in the text (section 4.1). 



in A(" +1 ), and the integrations are over those bonds b G dpW that also belong to the interior of 
The plaquette p is one of the tiling plaquettes (which one is immaterial by symmetry and 
translation invariance). Note that lattices A^ and A^ !+1 ' have lattice spacings differing by a scale 
factor of 2. To work out (4.6) explicitly, all factors are expanded in group characters, viz. 



exp(-S^[U[p^}}) 



c (n) +^0 j(n) X (U\p^}) 



exp 



A>(»)+E/Wz(tf[p W ]) 

j 



(4.7) 
(4.8) 



and similarly for [t/[p( n+1 )]] and the determinant factor. Inserting in (4.6), on the r.h.s. the 
various factor powers may be evaluated by the CG relations for combining characters and the inte- 
grations may be carried out exactly by using the character orthogonality relations. As a result one 
obtains explicit, albeit lengthy and ugly, recursion relations determining the expansion coefficients 
Cj(n + 1) at the (n + l)-th blocking step in terms of those at the n-the step. The procedure may 
then be iterated any number of times. It is important to note that the action at each step in general 
contains an infinite number of non-trivial characters, with corresponding couplings j3/(«) (cf. (4.8). 
They will be generic ally generated just after a single blocking step even if one starts with a single 
representation action, e.g the Wilson fundamental representation action, on the original lattice. 

This RG recursion scheme depends only on two quantities: the parameter T7, and, at each 
blocking step, the critical K(2"a). We have applied it to SU (2) and SU (3). The results are shown 
in Table 1 . Nt and Nu denote the number of flavors at which an IR fixed point first appears and 
that at which asymptotic freedom is lost, respectively. The parameter Tj was fixed by the conformal 
window upper bound value value A^u for SU (2) fundamental representation fermions. This deci- 
mation parameter could in principle be tuned independently in each case, but it turns out that once 
set by the top entry in the second column, no further adjustment was found necessary in order to 
produce the rest of the entries with the correct A\j values. Having fixed the parameter r\ the only 
uncertainties in practice are in obtaining the critical K values. They were obtained from estimates 
from (4.2) and/or values from simulations reported in the literature [14]. The resulting recursion 
relations for the character coefficients can then be run, very cheaply, to essentially arbitrary ac- 
curacy. It is noteworthy that such an approximate decimation scheme can already produce results 
such as shown in Table 1. 
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5. Discussion - Outlook 

Elucidating the phase diagram of gauge theories with different fermion content is a fully non- 
perturbative problem that presents many challenges. All methods that have been tried so far have 
their practical or conceptual shortcomings. Here we have considered what should be conceptu- 
ally the most straightforward approach, i.e. the flow of the system under successive RG blocking 
transformations. 

We have first explored locality-preserving RG blocking schemes for light fermions. We saw 
that blockings maintaining locality in an arbitrary gauge field background can be devised, but they 
generically require non-hypercubic blocked fermion actions. Such schemes may prove useful in 
investigating a variety of interesting issues and applications, e.g. the 'rooting' of fermion determi- 
nants; or improvements in the derivation of effective actions for exploring the QCD phase diagram 
via strong coupling and hopping expansions [15] - [16] combined with simulations. 

Full (fermion plus gauge field) RG blocking would seem to provide the most direct and reliable 
way for exploring the IR phase structure with varying fermion content. Carrying out the gauge field 
blocking explicitly in the presence of the blocked fermions, however, presents new challenges. In 
particular, implementation of MCRG procedures, which have been applied in the pure gauge theory 
case, is no longer straightforward for the reasons outlined above. Much work remains to be done 
to sort out these issues. 

At the same time judiciously chosen simplifications lead to approximate schemes which prove, 
as we saw, surprisingly effective. We found, in particular, that combining fermion blocking with 
gauge blocking within a potential moving approximation yields very reasonable results compared 
to other computationally much more expensive methods that have being tried. It is interesting to 
consider why this might be so. As it is well-known in the pure gauge case potential moving gives 
rather accurate beta functions at weak coupling (within 2% for SU (2)). When fermions are coupled 
to the system, appropriate weighing of the fermion contribution by a decimation parameter (such 
as the parameter r\ above) apparently can preserve this feature. This ensures that the upper limits 
of conformal windows as predicted by perturbation theory are correctly reproduced. It generally 
sets the flow 'on the right path'. It is very important in this connection that this flow moves off 
into an infinite dimensional subspace of couplings (cf. (4.8)) - we know from MCRG studies [12], 
that including enough characters in the effective action is important for good behavior over a wider 
range of scales. Moving in a sufficiently large, and preferably infinite dimensional, coupling space 
is crucial for evading lattice artifacts such as critical points and boundaries that can appear when 
projecting into a small finite-dimensional coupling constant subspace. 

If one is interested only in locating fixed points, rather than in detailed comparisons of the be- 
havior of observables, following the RG flow provides a 'robust' method in the following sense. As 
long as the flow manages to be within the right domain of attraction, it will lead to the correspond- 
ing fixed point; high accuracy in obtaining an effective action flow is not required. Approximate 
blocking schemes that incorporate the right features necessary to achieve this are then perfectly 
adequate for the purpose. This appears to indeed be the case for the scheme considered here. This 
minimal scheme, and other similarly constructed decimation schemes, are amenable to refinement 
by the inclusion of additional decimation parameters in more elaborate blocking definitions. It 
remains to be seen what improvements, if any, can be achieved by such refinements. 
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